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Classical novae are important contributors to the abundances of key isotopes, sucli as tiie ra- 
dioactive ^*F, whose observation by satellite missions could provide constraints on nucleosynthesis 
models in novae. The ^^0(p,7)^*F reaction plays a critical role in the synthesis of both oxygen and 
fluorine isotopes but its reaction rate is not well determined because of the lack of experimental data 
at energies relevant to novae explosions. In this study, the reaction cross section has been measured 
directly for the flrst time in a wide energy range Ecm — 200 — 370 keV appropriate to hydrogen 
burning in classical novae. In addition, the _E=183 keV resonance strength, cJ7=1.67±0.12 ^eV, has 
been measured with the highest precision to date. The uncertainty on the ^^0(p,7)^*F reaction rate 
has been reduced by a factor of 4, thus leading to firmer constraints on accurate models of novae 
nucleosynthesis. 

PACS numbers: 26.20.Cd; 26.30.-k; 26.50.+X 



Classical novae, a frequent phenomenon in our Galaxy, 
are explained as thermonuclear explosions on the sur- 
face of white dwarf stars accreting hydrogen-rich mate- 
rial from less evolved companions in binary star systems 
[l| and have been proposed as a key source of ^"^C, ^^N, 
-'^''d^O, and i8,i9p isotopes in the Universe 0- In particu- 
lar, the short-lived radioisotope ^^F {ti/2 = 110 min) may 
provide a signature of novae outbursts through the detec- 
tion of 511 keV 7-ray emission from positron-electron an- 
nihilation following its /?"*" decay. Indeed, the observation 
of these 7 rays by satellite missions could put constraints 
on current nova models Hydrogen burning of ^^O is 
believed to play a key role on the destruction of ^''O and 
on the formation of ^®F, mainly through the compet- 
ing reactions i^0(p,7)i^F and i'^0(p,Q;)i^N. Thus, the 
thermonuclear rates of both reactions should be deter- 
mined with a high degree of accuracy directly in the en- 
ergy region of hydrogen burning in classical novae. In 
this Letter, we report on a four-fold improvement in the 
^''0(p,7)^^F reaction rate determination. We have mea- 



sured the "'^''0(p,7)"'^^F reaction cross section down to the 
lowest energies to date and within the Gamow window 
3 for peak temperatures T = 0.1 - 0.4 GK . 

At nova temperatures, the ^''0(p,7)^^F reaction rate is 
dominated by a direct-capture (DC) reaction mechanism 
despite the presence of two narrow resonances at E = 66 
and 183 keV above the proton threshold in ^*F [^, 0| 
(all energies are in the center-of-mass system, unless oth- 
erwise stated). In addition, non-resonant contributions 
arise from the low-energy tails of two broad resonances 
a,t E = 556.7 and 676.7 keV. A reliable determination 
of the total ^^0(p,7)^^F reaction rate thus requires the 
accurate knowledge of the individual energy dependence 
of both resonant and non-resonant contributions. 

Early investigations by Rolfs [3l reported a constant 
DC contribution to the S-factor (Sdc) 0, Hi in agree- 
ment with the four lowest data points measured at en- 
ergies E — 280 — 425 keV by prompt 7-ray detection 
0. It was later questioned by Fox et al. Q whether 
these data points were dominated by the DC process or 
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FIG. 1. Sample spectrum from prompt 7-ray measurements of the ^^0(p,7)^®F reaction at the E — 183 keV resonance. 
Primary transitions are labelled as R^_Ex, while secondary transitions are labelled by their energies. Some of the most 
prominent background transitions (dashed arrows) are also shown for comparison, together with a time-normalised background 
spectrum (grey area). 



affected by thepresence of the two broad resonance tails. 
Thus, in Ref. [5|, the Sdc was calculated using measured 
spectroscopic factors and a realistic Woods-Saxon poten- 
tial and found to be up to a factor of '^2.5 lower than 
that reported in Ref. [7|]. A subsequent measurement 
at ~ 257 — 470 keV led to an Sdc in agreement 
with the predictions by Fox et al. [3| and still about a 
factor of 2 lower than in Ref. Q- More recently, the to- 
tal (i.e., DC plus broad- resonance contributions) S factor 
was measured at = 260 - 470 keV using the DRAGON 
recoil separator at TRIUMF and found to be in fairly 
good agreement with values by Ref. 0] , but consistently 
higher than the total S factor reported in Ref. 5]. As 
for the E = 183 keV resonance strength, only two val- 
ues exist to date: ujj — (1.2 ± 0.2) x 10~^ eV, as de- 
termined by a prompt 7-ray measurement 0, and 
(jj"f = (2.2 ± 0.4) X 10~^ eV, as determined by the ac- 
tivation technique [ll[, disagree at the 95% confidence 
level. The origin of this discrepancy is not understood at 
present, but may be due in part to unobserved gamma 
transitions in Refs. 

[i,|l^ and/or an inappropriate sub- 
traction of the DC component in either Refs. or 



Refs. [ll|. Thus, the lack of experimental data at low 



energies and the largely unconstrained Snc-factor have 
so far precluded the accurate determination of the ther- 
monuclear rate for this important reaction. 

Here, we report on the results of a new and improved 
investigation of the ^^0(p,7)^^F reaction using both acti- 



vation and prompt 7-ray detection techniques to address 
a key source of discrepancy in existing results. Measure- 
ments were carried out at the 400kV LUNA accelera- 
tor [l3| of the Laboratory for Underground Nuclear As- 
trophysics (LUNA) facility, which offers significant im- 
provements in sensitivity thanks to its low-background 
environment [11]. A proton beam, with currents up to 
200 fiA on target, entered the target chamber through a 
liquid-nitrogen cooled copper-pipe biased to -300 V for 
secondary electrons suppression. The target was directly 
water cooled with de-ionised water. Targets were pre- 
pared by anodization of Ta backings (0.3 mm thick disks) 
in isotopically enriched water (66% in ^^O and 4% in 
^^0). Full details on target preparation and characteri- 
sation have been reported in lj|. The target thickness 
was closely monitored for signs of degradation during in- 
tense proton-beam bombardment by regularly measuring 
the thick-target yield profile of the narrow isolated res- 



onance at = 143 keV in i®0(p,7)i^F [ij. Prompt 7 
rays from the i^0(p,7)i**F reaction (Q=5606.5±0.5 keV) 
were detected using a large volume (115% relative effi- 
ciency) high-purity germanium (HPGe) detector placed 
at 1.5 cm from the target and surrounded by 5 cm lead. 
Both detector and target were at an angle of 55° with re- 
spect to the beam axis, with the detector's front face par- 
allel to the target surface. Energy calibration, full-energy 
peak efficiency, and total efficiency were determined tak- 
ing into account corrections for true-coincidence sum- 
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FIG. 2. Total astrophysical S factor as a function of center- 
of-mass energy for the ^^0(p,7)^*F reaction. Filled symbols 
refer to the activation (circles) and the prompt 7-ray (squares) 
measurements in the present study. The solid line is the best 
fit to our data, while the dotted line labeled "Newton" is 
calculated using the DC contribution from Ref. and the 
resonance parameters from Ref. as used in Ref. [20l | (see 
text for further details). 



ming as described in Ref. [15|. Further details on the 
data analysis not reported^ here will be presented in a 
forthcoming publication 
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The ^''0(p,7)^*F reaction proceeds by populating sev- 
eral states in ^^F, leading to a complex decay scheme. 
Measurements were taken at the 183 keV resonance and 
at several ofF-resonance energies. At LUNA it was possi- 
ble to observe and identify a large number of the 7-ray 
transitions from the entrance channel to ^®F intermediate 
or ground state (the primaries) and subsequent direct or 
multistep decay of excited states to the "'^^F ground state 
(the secondaries) . A sample 7-ray spectrum showing the 
quality of our data is given in Fig. [TJ Independent deter- 
minations of the total S-factor were carried out for all pri- 
mary and secondary transitions. The main primary tran- 
sitions to i^F states at = 937,1121,3061,3791,3839 
and 4116 keV were analysed using the peak-shape ap- 
proach described in Ref. 17|. The resulting total astro- 



physical S-factor, obtained from the sum of all primary 
contributions, is shown in Fig. [H In total, 14 secondaries 
were also observed and analysed leading to an S factor 
(not shown in the Figure) in agreement with that from 
the primaries. Summing corrections amounted on aver- 
age to about 30% and were properly taken into account. 
In addition to the statistical errors shown, the primary 
total S-factor data carry a systematic uncertainty due 
to summing effects (1.5%), relative (1.5%) and absolute 
(3.5%) efficiency uncertainties, and an additional con- 
tribution of 7.6% (in common with the activation data) 
obtained as a sum in quadrature of uncertainties on col- 
lected charge (3%), stopping power (4%), stoichiometry 



(3%), target thickness (3.8%), and ^'^O abundance (3%). 

The £^=183 keV resonance strength ujj was deter- 
mined using the thick-target yield approach [4;]. Prior 
to our investigation the £^=183 keV resonance had been 
observed to decay to only two states at i<^x=1080 and 
937 keV in with branching ratios of 60% and 40%, 
respectively [18|. However, several additional transi- 
tions were identified in our study, as shown in Fig. [T] 
Both the resonance strength and the branching ratios 
of all observed transitions were treated as free param- 
eters in a g lobal St to the experimental counts as 



in Ref. [17[. For the three most dominant transitions. 



R^1080, R^-937 and R-^2101 keV, the branching ratios 
were determined as 41±3%, 25±2% and 12±2%, respec- 
tively. The branching ratios for a remaining six observed 
transitions were less than 6% each. Non-resonant con- 
tributions were subtracted for each transition individu- 
ally. The measured E=183 keV resonance strength was 
a;7=1.70±3%(stat)±7.7%(sys)AfeV. (Note that for mea- 
surements at the resonance, the target thickness uncer- 
tainty is only relevant for the DC subtraction and the 
common uncertainties reduce to 6.6%). 

An independent determination of both resonant and 
non-resonant contributions to the ^^0(p,7)^®F total cross 
section was also obtained with the activation technique, 
i.e. by detecting the 511 keV 7 rays from the positron 
annihilation following the decay of ^^F. Targets were 
irradiated for several hours to saturate the ^^F activ- 
ity, using the same setup described earlier but without 
lead shielding. Loss of ^^F from the target during irra- 
diation was estimated to be below 1% at all bombarding 
energies according to GEANT4 Monte Carlo simulations. 
Activation spectra were recorded every 20 minutes over a 
period of several hours using the low-background facility 
STELLA (SubTErranean Low Level Assay) |19|]. The ab- 
solute efficiency of the detector was measured with a cal- 
ibrated ®^Sr source which emits a single 7-ray line at 514 
keV allowing for a direct determination of the efficiency 
at the required energy, free from summing corrections. 
The laboratory background was negligible compared to 
the count rate of the activated targets. The irradiated 
targets were placed on top of the detector crystal and a 
Ta absorber, in direct contact with the activated target, 
was used to fully stop the emitted positrons. Our best 



value of t 



1/2 



109.6±0.4 min for the ^^F half-hfe was ob- 



tained as an average of several fits to experimental data 
at different energies and found to be in excellent agree- 
ment with the literature value ti/2 = 109.77 ± 0.05 min. 
The only contaminant to the spectra was found to be 
the short-lived ^"^N emitter. No longer-lived positron 
emitters were observed after several half-lives. From 
the number of observed disintegrations at off-resonance 
energies, the total S factor was determined following 
a standard procedure described in Ref. Results 
from the activation analysis are shown in Fig. [5] In 
addition to the statistical errors shown, data are af- 
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TABLE I. Reaction rate for the ^^0(p,7)^*I 
function of temperature. 

Temperature lower limit recommended value upper limit 
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fected by an independent systematic uncertainty of 2.9% 
(sum in quadrature of uncertainties from backscattering 
losses (1%), ^"^N contamination (1%), detector efficiency 
(2.5%)), and by a 7.6% uncertainty in common with tlie 
prompt 7-ray data (see above). Tlie E=183 keV reso- 
nance strength was derived from measurements at ener- 
gies on the plateau of the thick-target yield after subtrac- 
tion of a DC contribution of about 18%, which results in a 
3% DC-correction uncertainty on ujj. The obtained res- 
onance strength, a;7=1.65±1.8%(stat)±7.8%(sys) ^eV, 
is in excellent agreement with the result of the prompt 
7-ray measurements. The weighted average of the res- 
onance strength values from both measurements gives 
a;7=1.67±0.12(stat-|-sys) /ieV. 

To obtain an overall total S factor, the activation and 
the prompt 7-ray data sets were analyzed in a common 
fit procedure. Unfortunately, inclusion of other data sets 
(e.g. from Ref 0) in our fit was not possible since 
these seem to be affected by a mixture of systematic 
and point-to-point uncertainties that prevent a common 
and self-consistent analysis of all data sets. Following 
the phenomenological approach of Refs. , the fit- 

ting function included contributions from the two broad 
resonances at E=557 and 667 keV, with energy depen- 
dence of the involved partial widths described by the 
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FIG. 3. (Colour online) ^^0(p,7)^*F reaction rate ratio as 
a function of temperature. The (red) solid line is the ratio 
between the rate from a recent compilation rate (Ref. [20I ]) 
and the present value. Hatched and shaded areas represent a 
1(7 uncertainty on present and previous rate, respectively. In 
the relevant nova temperatures T = 0.1 — 0.4 GK, the rate 
uncertainty is reduced by a factor of 4 in the present study. 



broad-resonance formalism Q, as well as a constant DC 
term. The weak energy dependence of the DC compo- 
nent was neglected in the present analysis, as in previ- 
ous works [5,, 201 , since its effect is beyond the precision 
of our data. Finally, preliminary calculations 16[ indi- 
cate that interference effects between resonant and DC 
contributions are negligible and have not been included. 
A proper treatment of the systematic uncertainties 21 1 
required the introduction of three scaling factors, one 
for each data set to account for non-common system- 
atic uncertainties, and one for the common uncertain- 
ties. Their values were then obtained simultaneously in 
a modified fit (for a similar approach see Ref. 22]). 
Additional free parameters in the fit were the DC con- 
tribution and the partial F^ width of the 557 keV res- 
onance. The former parameter was left unconstrained 
while the latter was treated as a free parameter weighed 
by a Gaussian probability distribution with expectation 
value and standard deviation given by the literature value 
r-,,557 = 0.57 ± 0.13 eV [13]. The influence of the second 
broad-resonance on the fit quality was extremely small 
and, thus, F^.677 was fixed to the literature value f2^. 
The best fit (x^ = 13.9 for 19 data points and 5 fit pa- 
rameters) was obtained for F-y^557 = 0.70 eV and a DC 
component of Sdc = 4.4±0.4 keVb with scaling factors of 
Cact = 1.008, Cpiim — 0.972, and Ccom = 0.963. These fac- 
tors are in good agreement with their expectation value 
of 1, within the corresponding systematic uncertainties 
they represent. 

The astrophysical ^^0(p,7)^^F reaction rate was cal- 
culated using the formalism described in Ref. [3] as an 
incoherent sum of contributions from two narrow reso- 
nances (the E = 183 keV one from the present work, 
the E = 66 keV one from the literature [2^) and the 
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combined contribution of the broad resonances and the 
direct capture, this latter obtained as numerical integra- 
tion of the present S-factor best fit (solid curve in Fig. 

. The resulting reaction rate is tabulated in Table U 
with lower and upper limits given by the 68% confidence 
level. A comparison between our rate and that from a 
recent compilation [2^ (Fig. [S]) shows an improvement 
of a factor of 4 in the present rate uncertainty. 

On the basis of our improved rate, we have explored 
preliminary implications on the abundances of key iso- 
topes produced by classical novae. For example, ac- 
cording to Ref. [2^ a 25% (la) uncertainty in the 
^''0(p,7)^^F reaction rate leads to a 40 — 50% varia- 
tion on the calculated yields of ^^O, ^*F and ^^F at 
T = 0.1 - 0.4 GK. With the present rate being deter- 
mined at the 5% level (see Fig. [3]), the uncertainty on 
the ^^O, ^^F and ^^F abundances reduces to less than 
10% and puts firmer constraints towards more accurate 
nucleosynthesis calculations in novae events. The effects 
of our revised rate on the computation of detailed nova 
models will be discussed in a forthcoming paper. 

In summary, the ^^0(p,7)^*F reaction cross section has 
been measured for the first time in the relevant energy 
region of hydrogen burning in classical novae. The astro- 
physical S factor has been determined both by means of 
prompt 7-ray and activation measurements and results 
from the two approaches were found to be in excellent 
agreement. In addition, we have obtained the most accu- 
rate determination of the £^=183 keV resonance strength 
to date as W7=1.67±0.12 /leV. The reaction rate uncer- 
tainty was reduced by a factor of 4 leading to improved 
constraints on classical novae nucleosynthesis. 
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